home *** CD-ROM | disk | FTP | other *** search
| Unknown | 1991-01-05 | 9.8 KB |
open in: 
MacOS 8.1
|
Win98
|
DOS
view JSON data
|
view as text
This file was not able to be converted.
This format is not currently supported by dexvert.
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| file
| data
| default
| |
| 100%
| xdgMime
| image/x-tga
| default
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 0d 00 0a 15 f4 20 50 72 | 6f 67 72 61 6d 20 5a 65 |..... Pr|ogram Ze|
|00000010| 72 46 75 6e 63 0d 00 14 | 20 f4 20 28 43 29 20 44 |rFunc...| . (C) D|
|00000020| 2e 45 2e 20 42 6f 77 65 | 72 20 20 34 20 4a 61 6e |.E. Bowe|r 4 Jan|
|00000030| 20 31 39 39 31 0d 00 1e | 28 f4 20 56 65 72 73 69 | 1991...|(. Versi|
|00000040| 6f 6e 20 32 20 66 6f 72 | 20 41 72 63 68 69 6d 65 |on 2 for| Archime|
|00000050| 64 65 73 20 2b 20 42 41 | 53 49 43 20 56 0d 00 28 |des + BA|SIC V..(|
|00000060| 1e f4 20 28 63 6f 6e 76 | 65 72 74 65 64 20 66 72 |.. (conv|erted fr|
|00000070| 6f 6d 20 42 41 53 49 43 | 20 32 29 0d 00 32 05 20 |om BASIC| 2)..2. |
|00000080| 0d 00 3c 23 f4 20 52 65 | 61 6c 20 61 6e 64 2f 6f |..<#. Re|al and/o|
|00000090| 72 20 63 6f 6d 70 6c 65 | 78 20 7a 65 72 6f 65 73 |r comple|x zeroes|
|000000a0| 20 6f 66 0d 00 46 22 f4 | 20 61 6e 20 61 72 62 69 | of..F".| an arbi|
|000000b0| 74 72 61 72 79 20 66 75 | 6e 63 74 69 6f 6e 20 6f |trary fu|nction o|
|000000c0| 66 20 6f 6e 65 0d 00 50 | 23 f4 20 76 61 72 69 61 |f one..P|#. varia|
|000000d0| 62 6c 65 20 62 79 20 4d | 75 6c 6c 65 72 27 73 20 |ble by M|uller's |
|000000e0| 6d 65 74 68 6f 64 20 2e | 0d 00 5a 05 20 0d 00 64 |method .|..Z. ..d|
|000000f0| 3b f4 20 44 2e 45 2e 20 | 4d 75 6c 6c 65 72 20 22 |;. D.E. |Muller "|
|00000100| 41 20 6d 65 74 68 6f 64 | 20 6f 66 20 73 6f 6c 76 |A method| of solv|
|00000110| 69 6e 67 20 61 6c 67 65 | 62 72 61 69 63 20 65 71 |ing alge|braic eq|
|00000120| 75 61 74 69 6f 6e 73 22 | 0d 00 6e 24 f4 20 4d 54 |uations"|..n$. MT|
|00000130| 41 43 20 56 6f 6c 2e 20 | 31 30 20 70 70 20 32 30 |AC Vol. |10 pp 20|
|00000140| 38 2d 32 31 35 20 28 31 | 39 35 36 29 0d 00 78 04 |8-215 (1|956)..x.|
|00000150| 0d 00 82 35 f4 20 43 6f | 6e 74 65 20 26 20 64 65 |...5. Co|nte & de|
|00000160| 20 42 6f 6f 72 20 22 45 | 6c 65 6d 65 6e 74 61 72 | Boor "E|lementar|
|00000170| 79 20 4e 75 6d 65 72 69 | 63 61 6c 20 41 6e 61 6c |y Numeri|cal Anal|
|00000180| 79 73 69 73 22 0d 00 8c | 22 f4 20 32 6e 64 20 65 |ysis"...|". 2nd e|
|00000190| 64 69 74 69 6f 6e 20 4d | 63 47 72 61 77 2d 48 69 |dition M|cGraw-Hi|
|000001a0| 6c 6c 20 31 39 37 32 0d | 00 96 05 20 0d 00 a0 04 |ll 1972.|... ....|
|000001b0| 0d 00 aa 26 ee 20 85 20 | 3a 20 ee 20 85 20 87 20 |...&. . |: . . . |
|000001c0| 3a 20 40 25 3d 31 30 20 | 3a 20 ef 32 30 20 3a 20 |: @%=10 |: .20 : |
|000001d0| e5 20 8d 54 62 41 0d 00 | b4 05 20 0d 00 be 07 eb |. .TbA..|.. .....|
|000001e0| 20 37 0d 00 c8 09 f2 69 | 6e 66 6f 0d 00 d2 07 eb | 7.....i|nfo.....|
|000001f0| 20 31 0d 00 dc 0f f2 69 | 6e 69 74 69 61 6c 69 73 | 1.....i|nitialis|
|00000200| 65 0d 00 e6 0d f2 76 61 | 6c 69 64 61 74 65 0d 00 |e.....va|lidate..|
|00000210| f0 0d e7 20 76 61 6c 69 | 64 20 8c 0d 00 fa 10 f2 |... vali|d ......|
|00000220| 73 63 72 65 65 6e 5f 69 | 6e 69 74 0d 01 04 0b f2 |screen_i|nit.....|
|00000230| 73 65 61 72 63 68 0d 01 | 0e 05 cd 0d 01 18 04 0d |search..|........|
|00000240| 01 22 19 eb 20 30 20 3a | 20 fb 20 30 20 3a 20 fb |.".. 0 :| . 0 : .|
|00000250| 20 31 32 39 20 3a 20 db | 0d 01 2c 04 0d 01 36 10 | 129 : .|..,...6.|
|00000260| e7 20 ac 28 76 61 6c 69 | 64 29 20 8c 0d 01 40 26 |. .(vali|d) ...@&|
|00000270| f1 20 27 27 22 20 20 49 | 6e 70 75 74 20 44 61 74 |. ''" I|nput Dat|
|00000280| 61 20 69 73 20 4e 4f 54 | 20 56 61 6c 69 64 20 21 |a is NOT| Valid !|
|00000290| 21 22 0d 01 4a 2a f1 20 | 27 27 22 20 20 50 6c 65 |!"..J*. |''" Ple|
|000002a0| 61 73 65 20 43 68 65 63 | 6b 20 61 6e 64 20 43 6f |ase Chec|k and Co|
|000002b0| 72 72 65 63 74 20 44 61 | 74 61 22 27 0d 01 54 05 |rrect Da|ta"'..T.|
|000002c0| cd 0d 01 5e 27 e7 20 9f | 20 80 20 9f 3c 3e 31 37 |...^'. .| . .<>17|
|000002d0| 20 80 20 9e 20 8c 20 f6 | 3a f1 20 22 20 61 74 20 | . . . .|:. " at |
|000002e0| 6c 69 6e 65 20 22 3b 9e | 0d 01 68 04 0d 01 72 05 |line ";.|..h...r.|
|000002f0| e0 0d 01 7c 04 0d 01 86 | 3e 2a 2a 2a 2a 2a 2a 2a |...|....|>*******|
|00000300| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000310| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000320| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000330| 2a 2a 2a 0d 01 90 04 0d | 01 9a 0b dd 20 f2 69 6e |***.....|.... .in|
|00000340| 66 6f 0d 01 a4 0b 76 61 | 6c 69 64 3d b9 0d 01 ae |fo....va|lid=....|
|00000350| 4b e3 20 72 6f 77 3d 30 | 20 b8 20 31 3a f1 8a 30 |K. row=0| . 1:..0|
|00000360| 2c 72 6f 77 29 3b bd 31 | 34 31 3b bd 31 33 30 3b |,row);.1|41;.130;|
|00000370| 22 43 6f 6d 70 6c 65 78 | 20 5a 65 72 6f 65 73 20 |"Complex| Zeroes |
|00000380| 6f 66 20 41 72 62 69 74 | 72 61 72 79 20 46 75 6e |of Arbit|rary Fun|
|00000390| 63 74 69 6f 6e 22 3a ed | 0d 01 b8 2e f1 27 22 20 |ction":.|.....'" |
|000003a0| 20 50 72 6f 67 72 61 6d | 20 5a 65 72 46 75 6e 63 | Program| ZerFunc|
|000003b0| 20 73 65 61 72 63 68 65 | 73 20 66 6f 72 20 61 6c | searche|s for al|
|000003c0| 6c 20 74 68 65 22 0d 01 | c2 2d f1 20 22 20 20 72 |l the"..|.-. " r|
|000003d0| 65 61 6c 20 6f 72 20 63 | 6f 6d 70 6c 65 78 20 7a |eal or c|omplex z|
|000003e0| 65 72 6f 65 73 20 6f 66 | 20 61 20 67 65 6e 65 72 |eroes of| a gener|
|000003f0| 61 6c 22 0d 01 cc 2d f1 | 20 22 20 20 66 75 6e 63 |al"...-.| " func|
|00000400| 74 69 6f 6e 20 6f 66 20 | 74 68 65 20 63 6f 6d 70 |tion of |the comp|
|00000410| 6c 65 78 20 76 61 72 69 | 61 62 6c 65 20 5a 2e 22 |lex vari|able Z."|
|00000420| 0d 01 d6 2c f1 27 22 20 | 20 49 74 20 75 73 65 73 |...,.'" | It uses|
|00000430| 20 74 68 65 20 4d 75 6c | 6c 65 72 20 61 6c 67 6f | the Mul|ler algo|
|00000440| 72 69 74 68 6d 20 77 68 | 69 63 68 22 0d 01 e0 2e |rithm wh|ich"....|
|00000450| f1 20 22 20 20 72 65 70 | 65 61 74 65 64 6c 79 20 |. " rep|eatedly |
|00000460| 63 6f 6d 70 75 74 65 73 | 20 61 6e 20 61 70 70 72 |computes| an appr|
|00000470| 6f 78 69 6d 61 74 69 6f | 6e 22 0d 01 ea 2e f1 20 |oximatio|n"..... |
|00000480| 22 20 20 74 6f 20 74 68 | 65 20 6e 65 78 74 20 72 |" to th|e next r|
|00000490| 6f 6f 74 20 62 79 20 65 | 78 74 72 61 70 6f 6c 61 |oot by e|xtrapola|
|000004a0| 74 69 6f 6e 20 6f 66 22 | 0d 01 f4 2e f1 20 22 20 |tion of"|..... " |
|000004b0| 20 61 20 70 61 72 61 62 | 6f 6c 61 20 66 69 74 20 | a parab|ola fit |
|000004c0| 74 6f 20 74 68 65 20 74 | 68 72 65 65 20 70 72 65 |to the t|hree pre|
|000004d0| 76 69 6f 75 73 22 0d 01 | fe 2c f1 20 22 20 20 61 |vious"..|.,. " a|
|000004e0| 70 70 72 6f 78 69 6d 61 | 74 69 6f 6e 73 2e 20 54 |pproxima|tions. T|
|000004f0| 68 69 73 20 6d 65 74 68 | 6f 64 20 69 73 20 66 61 |his meth|od is fa|
|00000500| 72 22 0d 02 08 2d f1 20 | 22 20 20 6d 6f 72 65 20 |r"...-. |" more |
|00000510| 72 6f 62 75 73 74 20 74 | 68 61 6e 20 74 68 65 20 |robust t|han the |
|00000520| 4e 65 77 74 6f 6e 2d 52 | 61 70 68 73 6f 6e 22 0d |Newton-R|aphson".|
|00000530| 02 12 2e f1 20 22 20 20 | 70 72 6f 63 65 64 75 72 |.... " |procedur|
|00000540| 65 20 61 6e 64 20 61 6c | 73 6f 20 68 61 73 20 6d |e and al|so has m|
|00000550| 75 63 68 20 73 75 70 65 | 72 69 6f 72 22 0d 02 1c |uch supe|rior"...|
|00000560| 28 f1 20 22 20 20 67 6c | 6f 62 61 6c 20 63 6f 6e |(. " gl|obal con|
|00000570| 76 65 72 67 65 6e 63 65 | 20 70 72 6f 70 65 72 74 |vergence| propert|
|00000580| 69 65 73 2e 22 0d 02 26 | 2c f1 27 22 20 20 54 68 |ies."..&|,.'" Th|
|00000590| 69 73 20 64 65 6d 6f 20 | 69 73 20 73 65 74 20 75 |is demo |is set u|
|000005a0| 70 20 66 6f 72 20 74 68 | 65 20 7a 65 72 6f 65 73 |p for th|e zeroes|
|000005b0| 22 0d 02 30 2e f1 20 22 | 20 20 6f 66 20 74 68 65 |"..0.. "| of the|
|000005c0| 20 72 65 61 6c 20 70 6f | 6c 79 6e 6f 6d 69 61 6c | real po|lynomial|
|000005d0| 20 46 28 5a 29 3d 5a 5e | 33 30 20 2d 20 31 22 0d | F(Z)=Z^|30 - 1".|
|000005e0| 02 3a 2e f1 20 22 20 20 | 62 75 74 20 74 68 65 20 |.:.. " |but the |
|000005f0| 63 6f 6e 74 72 6f 6c 20 | 64 61 74 61 20 6d 61 79 |control |data may|
|00000600| 20 62 65 20 6d 6f 64 69 | 66 69 65 64 22 0d 02 44 | be modi|fied"..D|
|00000610| 30 f1 20 22 20 20 74 6f | 20 61 6e 61 6c 79 73 65 |0. " to| analyse|
|00000620| 20 61 20 75 73 65 72 2d | 73 75 70 70 6c 69 65 64 | a user-|supplied|
|00000630| 20 66 75 6e 63 74 69 6f | 6e 2e 22 27 27 0d 02 4e | functio|n."''..N|
|00000640| 2e f1 20 bd 31 33 32 3b | 22 20 20 20 20 20 20 20 |.. .132;|" |
|00000650| 20 20 50 72 65 73 73 20 | 61 20 6b 65 79 2e 2e 2e | Press |a key...|
|00000660| 2e 2e 22 3a 2a 46 58 31 | 35 2c 31 0d 02 58 07 47 |..":*FX1|5,1..X.G|
|00000670| 3d a5 0d 02 62 05 e1 0d | 02 6c 04 0d 02 76 04 0d |=...b...|.l...v..|
|00000680| 02 80 11 dd 20 f2 69 6e | 69 74 69 61 6c 69 73 65 |.... .in|itialise|
|00000690| 0d 02 8a 04 0d 02 94 12 | ef 20 31 39 2c 31 2c 34 |........|. 19,1,4|
|000006a0| 2c 30 2c 30 2c 30 0d 02 | 9e 44 ef 20 32 34 2c 30 |,0,0,0..|.D. 24,0|
|000006b0| 3b 30 3b 38 30 30 3b 31 | 30 32 33 3b 20 20 3a 20 |;0;800;1|023; : |
|000006c0| f4 20 55 73 65 20 6c 65 | 66 74 20 20 70 6f 72 74 |. Use le|ft port|
|000006d0| 69 6f 6e 20 6f 66 20 73 | 63 72 65 65 6e 20 66 6f |ion of s|creen fo|
|000006e0| 72 20 67 72 61 70 68 69 | 63 73 0d 02 a8 44 ef 20 |r graphi|cs...D. |
|000006f0| 32 38 2c 32 36 2c 33 31 | 2c 33 39 2c 30 20 20 20 |28,26,31|,39,0 |
|00000700| 20 20 3a 20 f4 20 55 73 | 65 20 72 69 67 68 74 20 | : . Us|e right |
|00000710| 70 6f 72 74 69 6f 6e 20 | 6f 66 20 73 63 72 65 65 |portion |of scree|
|00000720| 6e 20 66 6f 72 20 74 65 | 78 74 20 6f 2f 70 0d 02 |n for te|xt o/p..|
|00000730| b2 39 ef 20 32 33 3b 38 | 32 30 32 3b 30 3b 30 3b |.9. 23;8|202;0;0;|
|00000740| 30 3b 20 20 20 20 3a 20 | f4 20 41 6e 63 69 65 6e |0; : |. Ancien|
|00000750| 74 20 6d 65 74 68 6f 64 | 20 66 6f 72 20 63 75 72 |t method| for cur|
|00000760| 73 6f 72 20 6f 66 66 0d | 02 bc 05 20 0d 02 c6 3d |sor off.|... ...=|
|00000770| 64 65 62 75 67 3d 30 20 | 20 3a 20 f4 20 53 65 74 |debug=0 | : . Set|
|00000780| 20 74 6f 20 31 20 66 6f | 72 20 70 72 69 6e 74 65 | to 1 fo|r printe|
|00000790| 72 20 6f 2f 70 20 6f 72 | 20 32 20 66 6f 72 20 73 |r o/p or| 2 for s|
|000007a0| 63 72 65 65 6e 20 6f 2f | 70 0d 02 d0 22 45 31 3d |creen o/|p..."E1=|
|000007b0| 31 45 2d 38 20 20 3a 20 | f4 20 52 6f 6f 74 20 50 |1E-8 : |. Root P|
|000007c0| 6f 73 69 74 69 6f 6e 20 | 54 6f 6c 0d 02 da 1f 45 |osition |Tol....E|
|000007d0| 32 3d 31 45 2d 38 20 20 | 3a 20 f4 20 46 75 6e 63 |2=1E-8 |: . Func|
|000007e0| 20 56 61 6c 75 65 20 54 | 6f 6c 0d 02 e4 22 45 33 | Value T|ol..."E3|
|000007f0| 3d 31 45 2d 31 34 20 3a | 20 f4 20 4d 75 6c 74 69 |=1E-14 :| . Multi|
|00000800| 70 6c 65 20 72 6f 6f 74 | 20 54 6f 6c 0d 02 ee 20 |ple root| Tol... |
|00000810| 45 34 3d 31 45 2d 31 36 | 20 3a 20 f4 20 43 53 51 |E4=1E-16| : . CSQ|
|00000820| 52 54 20 41 62 6f 72 74 | 20 54 6f 6c 0d 02 f8 25 |RT Abort| Tol...%|
|00000830| 45 31 39 3d 31 45 31 39 | 20 3a 20 f4 20 43 53 51 |E19=1E19| : . CSQ|
|00000840| 52 54 20 6f 76 65 72 66 | 6c 6f 77 20 6c 69 6d 69 |RT overf|low limi|
|00000850| 74 0d 03 02 0a 45 36 3d | 31 45 36 0d 03 0c 32 65 |t....E6=|1E6...2e|
|00000860| 70 73 3d 31 45 2d 33 20 | 3a 20 f4 20 50 65 72 74 |ps=1E-3 |: . Pert|
|00000870| 75 72 62 61 74 69 6f 6e | 20 66 6f 72 20 69 6e 69 |urbation| for ini|
|00000880| 74 69 61 6c 20 65 73 74 | 69 6d 61 74 65 0d 03 16 |tial est|imate...|
|00000890| 20 4d 61 78 69 74 3d 31 | 30 30 3a 20 f4 20 4d 61 | Maxit=1|00: . Ma|
|000008a0| 78 20 69 74 65 72 20 2f | 20 72 6f 6f 74 0d 03 20 |x iter /| root.. |
|000008b0| 1b 42 6c 61 6e 6b 24 3d | 22 20 20 20 20 20 20 20 |.Blank$=|" |
|000008c0| 20 20 20 20 20 20 20 22 | 0d 03 2a 04 0d 03 34 05 | "|..*...4.|
|000008d0| e1 0d 03 3e 04 0d 03 48 | 04 0d 03 52 0f dd 20 f2 |...>...H|...R.. .|
|000008e0| 76 61 6c 69 64 61 74 65 | 0d 03 5c 04 0d 03 66 17 |validate|..\...f.|
|000008f0| f3 20 46 24 2c 4e 72 2c | 43 24 20 3a 20 66 6e 24 |. F$,Nr,|C$ : fn$|
|00000900| 3d 46 24 0d 03 70 27 e7 | 20 46 24 3c 3e 22 50 4f |=F$..p'.| F$<>"PO|
|00000910| 4c 59 22 20 80 20 46 24 | 3c 3e 22 46 55 4e 43 22 |LY" . F$|<>"FUNC"|
|00000920| 20 8c 20 76 61 6c 69 64 | 3d a3 0d 03 7a 0d 46 24 | . valid|=...z.F$|
|00000930| 3d c0 46 24 2c 31 29 0d | 03 84 20 e7 20 4e 72 3c |=.F$,1).|.. . Nr<|
|00000940| 3e a8 28 4e 72 29 20 84 | 20 4e 72 3c 31 20 8c 20 |>.(Nr) .| Nr<1 . |
|00000950| 76 61 6c 69 64 3d a3 0d | 03 8e 1b de 20 5a 72 31 |valid=..|.... Zr1|
|00000960| 28 4e 72 29 2c 5a 72 32 | 28 4e 72 29 2c 43 28 4e |(Nr),Zr2|(Nr),C(N|
|00000970| 72 29 0d 03 98 2a e7 20 | 43 24 3c 3e 22 52 45 41 |r)...*. |C$<>"REA|
|00000980| 4c 22 20 80 20 43 24 3c | 3e 22 43 4f 4d 50 4c 45 |L" . C$<|>"COMPLE|
|00000990| 58 22 20 8c 20 76 61 6c | 69 64 3d a3 0d 03 a2 21 |X" . val|id=....!|
|000009a0| e7 20 43 24 3d 22 52 45 | 41 4c 22 20 63 6d 70 6c |. C$="RE|AL" cmpl|
|000009b0| 78 3d 30 20 8b 20 63 6d | 70 6c 78 3d 31 0d 03 ac |x=0 . cm|plx=1...|
|000009c0| 19 f3 20 58 6d 69 6e 2c | 58 6d 61 78 2c 59 6d 69 |.. Xmin,|Xmax,Ymi|
|000009d0| 6e 2c 59 6d 61 78 0d 03 | b6 27 e7 20 58 6d 69 6e |n,Ymax..|.'. Xmin|
|000009e0| 3e 3d 58 6d 61 78 20 84 | 20 59 6d 69 6e 3e 3d 59 |>=Xmax .| Ymin>=Y|
|000009f0| 6d 61 78 20 8c 20 76 61 | 6c 69 64 3d a3 0d 03 c0 |max . va|lid=....|
|00000a00| 0e e7 20 46 24 3d 22 50 | 22 20 8c 0d 03 ca 34 f2 |.. F$="P|" ....4.|
|00000a10| 63 6f 65 66 66 20 20 20 | 20 20 20 20 20 20 20 3a |coeff | :|
|00000a20| 20 f4 20 72 65 61 64 20 | 70 6f 6c 79 6e 6f 6d 69 | . read |polynomi|
|00000a30| 61 6c 20 63 6f 65 66 66 | 69 63 69 65 6e 74 73 0d |al coeff|icients.|
|00000a40| 03 d4 43 4d 61 78 6d 6f | 64 3d 31 45 33 38 5e 28 |..CMaxmo|d=1E38^(|
|00000a50| 32 2f 4e 72 29 20 3a 20 | f4 20 65 73 74 69 6d 61 |2/Nr) : |. estima|
|00000a60| 74 65 20 6d 61 78 20 73 | 65 6e 73 69 62 6c 65 20 |te max s|ensible |
|00000a70| 69 6e 70 75 74 20 74 6f | 20 70 6f 6c 79 20 65 76 |input to| poly ev|
|00000a80| 61 6c 0d 03 de 05 cd 0d | 03 e8 04 0d 03 f2 05 e1 |al......|........|
|00000a90| 0d 03 fc 04 0d 04 06 04 | 0d 04 10 0d dd 20 f2 73 |........|..... .s|
|00000aa0| 65 61 72 63 68 0d 04 1a | 05 20 0d 04 24 11 e3 20 |earch...|. ..$.. |
|00000ab0| 72 6f 6f 74 3d 31 20 b8 | 20 4e 72 0d 04 2e 04 0d |root=1 .| Nr.....|
|00000ac0| 04 38 3a 52 31 3d 30 20 | 3a 20 52 32 3d 30 20 3a |.8:R1=0 |: R2=0 :|
|00000ad0| 20 49 74 65 72 3d 30 20 | 3a 20 f4 20 53 65 74 20 | Iter=0 |: . Set |
|00000ae0| 69 6e 69 74 69 61 6c 20 | 65 73 74 69 6d 61 74 65 |initial |estimate|
|00000af0| 20 74 6f 20 28 30 2c 30 | 29 0d 04 42 04 0d 04 4c | to (0,0|)..B...L|
|00000b00| 05 f5 0d 04 56 04 0d 04 | 60 17 48 31 3d 2e 34 20 |....V...|`.H1=.4 |
|00000b10| 3a 20 48 32 3d 2e 34 2a | 63 6d 70 6c 78 0d 04 6a |: H2=.4*|cmplx..j|
|00000b20| 23 4c 61 6d 64 61 31 3d | 2d 30 2e 35 20 3a 20 4c |#Lamda1=|-0.5 : L|
|00000b30| 61 6d 64 61 32 3d 30 20 | 3a 20 63 6f 6e 76 3d 30 |amda2=0 |: conv=0|
|00000b40| 0d 04 74 04 0d 04 7e 2e | f4 20 46 69 72 73 74 20 |..t...~.|. First |
|00000b50| 74 68 72 65 65 20 65 73 | 74 69 6d 61 74 65 73 20 |three es|timates |
|00000b60| 61 72 65 20 20 52 2b 48 | 20 2c 20 52 2d 48 20 2c |are R+H| , R-H ,|
|00000b70| 20 52 0d 04 88 04 0d 04 | 92 05 f5 0d 04 9c 04 0d | R......|........|
|00000b80| 04 a6 0c c8 8e 20 65 72 | 72 20 ca 0d 04 b0 12 c9 |..... er|r ......|
|00000b90| 20 31 3a 20 52 31 3d 52 | 31 2b 65 70 73 0d 04 ba | 1: R1=R|1+eps...|
|00000ba0| 20 c9 20 32 3a 20 52 31 | 3d 2e 32 2a b3 2d 2e 31 | . 2: R1|=.2*.-.1|
|00000bb0| 20 3a 20 52 32 3d 2e 32 | 2a b3 2d 2e 31 0d 04 c4 | : R2=.2|*.-.1...|
|00000bc0| 05 cb 0d 04 ce 04 0d 04 | d8 17 5a 31 3d 52 31 2b |........|..Z1=R1+|
|00000bd0| 48 31 20 3a 20 5a 32 3d | 52 32 2b 48 32 0d 04 e2 |H1 : Z2=|R2+H2...|
|00000be0| 09 f2 65 76 61 6c 0d 04 | ec 05 20 0d 04 f6 0d e7 |..eval..|.. .....|
|00000bf0| 20 65 72 72 3d 30 20 8c | 0d 05 00 13 46 6d 31 3d | err=0 .|....Fm1=|
|00000c00| 46 31 20 3a 20 46 6d 32 | 3d 46 32 0d 05 0a 17 5a |F1 : Fm2|=F2....Z|
|00000c10| 31 3d 52 31 2d 48 31 20 | 3a 20 5a 32 3d 52 32 2d |1=R1-H1 |: Z2=R2-|
|00000c20| 48 32 0d 05 14 09 f2 65 | 76 61 6c 0d 05 1e 05 cd |H2.....e|val.....|
|00000c30| 0d 05 28 05 20 0d 05 32 | 0d e7 20 65 72 72 3d 30 |..(. ..2|.. err=0|
|00000c40| 20 8c 0d 05 3c 17 4f 6c | 64 66 31 3d 46 31 20 3a | ...<.Ol|df1=F1 :|
|00000c50| 20 4f 6c 64 66 32 3d 46 | 32 0d 05 46 1b 46 6d 31 | Oldf2=F|2..F.Fm1|
|00000c60| 3d 46 31 2d 46 6d 31 20 | 3a 20 46 6d 32 3d 46 32 |=F1-Fm1 |: Fm2=F2|
|00000c70| 2d 46 6d 32 0d 05 50 11 | 5a 31 3d 52 31 20 3a 20 |-Fm2..P.|Z1=R1 : |
|00000c80| 5a 32 3d 52 32 0d 05 5a | 09 f2 65 76 61 6c 0d 05 |Z2=R2..Z|..eval..|
|00000c90| 64 05 cd 0d 05 6e 05 20 | 0d 05 78 0b fd 20 65 72 |d....n. |..x.. er|
|00000ca0| 72 3d 30 0d 05 82 04 0d | 05 8c 47 f4 20 50 61 72 |r=0.....|..G. Par|
|00000cb0| 61 62 6f 6c 69 63 20 45 | 78 74 72 61 70 6f 6c 61 |abolic E|xtrapola|
|00000cc0| 74 69 6f 6e 20 2f 20 55 | 70 64 61 74 65 20 69 74 |tion / U|pdate it|
|00000cd0| 65 72 61 74 69 6f 6e 20 | 28 20 63 6f 6d 70 6c 65 |eration |( comple|
|00000ce0| 78 20 61 72 69 74 68 6d | 65 74 69 63 20 29 0d 05 |x arithm|etic )..|
|00000cf0| 96 36 f4 20 44 65 6c 66 | 3d 46 2d 4f 6c 64 66 20 |.6. Delf|=F-Oldf |
|00000d00| 3a 20 44 65 6c 74 3d 46 | 6d 2a 4c 61 6d 64 61 20 |: Delt=F|m*Lamda |
|00000d10| 3a 20 4e 75 6d 3d 2d 32 | 2a 46 2a 28 31 2b 4c 61 |: Num=-2|*F*(1+La|
|00000d20| 6d 64 61 29 0d 05 a0 23 | f4 20 47 3d 28 31 2b 32 |mda)...#|. G=(1+2|
|00000d30| 2a 4c 61 6d 64 61 29 2a | 44 65 6c 66 2d 4c 61 6d |*Lamda)*|Delf-Lam|
|00000d40| 64 61 2a 44 65 6c 74 0d | 05 aa 31 f4 20 42 3d 47 |da*Delt.|..1. B=G|
|00000d50| 2a 47 20 2b 20 32 2a 4e | 75 6d 2a 4c 61 6d 64 61 |*G + 2*N|um*Lamda|
|00000d60| 2a 28 44 65 6c 66 2d 44 | 65 6c 74 29 20 3a 20 53 |*(Delf-D|elt) : S|
|00000d70| 71 3d 53 51 52 28 42 29 | 0d 05 b4 2c f4 20 44 65 |q=SQR(B)|...,. De|
|00000d80| 6e 3d 47 20 2b 2d 20 53 | 71 20 28 46 6f 72 63 65 |n=G +- S|q (Force|
|00000d90| 20 6c 61 72 67 65 72 20 | 64 65 6e 6f 6d 69 6e 61 | larger |denomina|
|00000da0| 74 6f 72 29 0d 05 be 27 | f4 20 4c 61 6d 64 61 3d |tor)...'|. Lamda=|
|00000db0| 4e 75 6d 2f 44 65 6e 20 | 3a 20 48 3d 48 2a 4c 61 |Num/Den |: H=H*La|
|00000dc0| 6d 64 61 20 3a 20 5a 3d | 5a 2b 48 0d 05 c8 04 0d |mda : Z=|Z+H.....|
|00000dd0| 05 d2 05 f5 0d 05 dc 05 | 20 0d 05 e6 0d 63 66 6c |........| ....cfl|
|00000de0| 67 3d 63 6f 6e 76 0d 05 | f0 0a 63 6f 6e 76 3d 30 |g=conv..|..conv=0|
|00000df0| 0d 05 fa 05 20 0d 06 04 | 12 44 65 6c 66 31 3d 46 |.... ...|.Delf1=F|
|00000e00| 31 2d 4f 6c 64 66 31 0d | 06 0e 12 44 65 6c 66 32 |1-Oldf1.|...Delf2|
|00000e10| 3d 46 32 2d 4f 6c 64 66 | 32 0d 06 18 21 44 65 6c |=F2-Oldf|2...!Del|
|00000e20| 74 31 3d 4c 61 6d 64 61 | 31 2a 46 6d 31 20 2d 20 |t1=Lamda|1*Fm1 - |
|00000e30| 4c 61 6d 64 61 32 2a 46 | 6d 32 0d 06 22 21 44 65 |Lamda2*F|m2.."!De|
|00000e40| 6c 74 32 3d 4c 61 6d 64 | 61 32 2a 46 6d 31 20 2b |lt2=Lamd|a2*Fm1 +|
|00000e50| 20 4c 61 6d 64 61 31 2a | 46 6d 32 0d 06 2c 0f 4c | Lamda1*|Fm2..,.L|
|00000e60| 31 3d 31 2b 4c 61 6d 64 | 61 31 0d 06 36 1f 4e 75 |1=1+Lamd|a1..6.Nu|
|00000e70| 6d 31 3d 2d 32 2a 28 4c | 31 2a 46 31 20 2d 20 4c |m1=-2*(L|1*F1 - L|
|00000e80| 61 6d 64 61 32 2a 46 32 | 29 0d 06 40 1f 4e 75 6d |amda2*F2|)..@.Num|
|00000e90| 32 3d 2d 32 2a 28 4c 31 | 2a 46 32 20 2b 20 4c 61 |2=-2*(L1|*F2 + La|
|00000ea0| 6d 64 61 32 2a 46 31 29 | 0d 06 4a 10 4c 31 3d 4c |mda2*F1)|..J.L1=L|
|00000eb0| 31 2b 4c 61 6d 64 61 31 | 0d 06 54 14 4c 32 3d 4c |1+Lamda1|..T.L2=L|
|00000ec0| 61 6d 64 61 32 2b 4c 61 | 6d 64 61 32 0d 06 5e 38 |amda2+La|mda2..^8|
|00000ed0| 47 31 3d 4c 31 2a 44 65 | 6c 66 31 20 2d 20 4c 32 |G1=L1*De|lf1 - L2|
|00000ee0| 2a 44 65 6c 66 32 20 2d | 20 4c 61 6d 64 61 31 2a |*Delf2 -| Lamda1*|
|00000ef0| 44 65 6c 74 31 20 2b 20 | 4c 61 6d 64 61 32 2a 44 |Delt1 + |Lamda2*D|
|00000f00| 65 6c 74 32 0d 06 68 38 | 47 32 3d 4c 32 2a 44 65 |elt2..h8|G2=L2*De|
|00000f10| 6c 66 31 20 2b 20 4c 31 | 2a 44 65 6c 66 32 20 2d |lf1 + L1|*Delf2 -|
|00000f20| 20 4c 61 6d 64 61 32 2a | 44 65 6c 74 31 20 2d 20 | Lamda2*|Delt1 - |
|00000f30| 4c 61 6d 64 61 31 2a 44 | 65 6c 74 32 0d 06 72 24 |Lamda1*D|elt2..r$|
|00000f40| 4e 31 3d 32 2a 28 4e 75 | 6d 31 2a 4c 61 6d 64 61 |N1=2*(Nu|m1*Lamda|
|00000f50| 31 20 2d 20 4e 75 6d 32 | 2a 4c 61 6d 64 61 32 29 |1 - Num2|*Lamda2)|
|00000f60| 0d 06 7c 24 4e 32 3d 32 | 2a 28 4e 75 6d 32 2a 4c |..|$N2=2|*(Num2*L|
|00000f70| 61 6d 64 61 31 20 2b 20 | 4e 75 6d 31 2a 4c 61 6d |amda1 + |Num1*Lam|
|00000f80| 64 61 32 29 0d 06 86 12 | 44 31 3d 44 65 6c 66 31 |da2)....|D1=Delf1|
|00000f90| 2d 44 65 6c 74 31 0d 06 | 90 12 44 32 3d 44 65 6c |-Delt1..|..D2=Del|
|00000fa0| 66 32 2d 44 65 6c 74 32 | 0d 06 9a 24 42 31 3d 47 |f2-Delt2|...$B1=G|
|00000fb0| 31 2a 47 31 20 2d 20 47 | 32 2a 47 32 20 2b 20 4e |1*G1 - G|2*G2 + N|
|00000fc0| 31 2a 44 31 20 2d 20 4e | 32 2a 44 32 0d 06 a4 24 |1*D1 - N|2*D2...$|
|00000fd0| 42 32 3d 20 20 20 20 32 | 2a 47 31 2a 47 32 20 20 |B2= 2|*G1*G2 |
|00000fe0| 20 2b 20 4e 32 2a 44 31 | 20 2b 20 4e 31 2a 44 32 | + N2*D1| + N1*D2|
|00000ff0| 0d 06 ae 05 20 0d 06 b8 | 09 f2 73 71 72 74 0d 06 |.... ...|..sqrt..|
|00001000| c2 05 20 0d 06 cc 1a 73 | 67 6e 3d b4 28 47 31 2a |.. ....s|gn=.(G1*|
|00001010| 53 71 31 20 2b 20 47 32 | 2a 53 71 32 29 0d 06 d6 |Sq1 + G2|*Sq2)...|
|00001020| 15 44 65 6e 31 3d 47 31 | 20 2b 20 73 67 6e 2a 53 |.Den1=G1| + sgn*S|
|00001030| 71 31 0d 06 e0 15 44 65 | 6e 32 3d 47 32 20 2b 20 |q1....De|n2=G2 + |
|00001040| 73 67 6e 2a 53 71 32 0d | 06 ea 1f 44 6e 6f 72 6d |sgn*Sq2.|...Dnorm|
|00001050| 3d 44 65 6e 31 2a 44 65 | 6e 31 20 2b 20 44 65 6e |=Den1*De|n1 + Den|
|00001060| 32 2a 44 65 6e 32 0d 06 | f4 23 e7 20 44 6e 6f 72 |2*Den2..|.#. Dnor|
|00001070| 6d 3d 30 20 44 65 6e 31 | 3d 31 3a 44 65 6e 32 3d |m=0 Den1|=1:Den2=|
|00001080| 30 3a 44 6e 6f 72 6d 3d | 31 0d 06 fe 04 0d 07 08 |0:Dnorm=|1.......|
|00001090| 16 f4 20 55 70 64 61 74 | 65 20 69 74 65 72 61 74 |.. Updat|e iterat|
|000010a0| 69 6f 6e 0d 07 12 04 0d | 07 1c 28 4c 61 6d 64 61 |ion.....|..(Lamda|
|000010b0| 31 3d 28 4e 75 6d 31 2a | 44 65 6e 31 20 2b 20 4e |1=(Num1*|Den1 + N|
|000010c0| 75 6d 32 2a 44 65 6e 32 | 29 2f 44 6e 6f 72 6d 0d |um2*Den2|)/Dnorm.|
|000010d0| 07 26 28 4c 61 6d 64 61 | 32 3d 28 4e 75 6d 32 2a |.&(Lamda|2=(Num2*|
|000010e0| 44 65 6e 31 20 2d 20 4e | 75 6d 31 2a 44 65 6e 32 |Den1 - N|um1*Den2|
|000010f0| 29 2f 44 6e 6f 72 6d 0d | 07 30 0c 4f 6c 64 66 31 |)/Dnorm.|.0.Oldf1|
|00001100| 3d 46 31 0d 07 3a 0c 4f | 6c 64 66 32 3d 46 32 0d |=F1..:.O|ldf2=F2.|
|00001110| 07 44 0d 46 6d 31 3d 44 | 65 6c 66 31 0d 07 4e 0d |.D.Fm1=D|elf1..N.|
|00001120| 46 6d 32 3d 44 65 6c 66 | 32 0d 07 58 1c 54 20 3d |Fm2=Delf|2..X.T =|
|00001130| 48 31 2a 4c 61 6d 64 61 | 31 20 2d 20 48 32 2a 4c |H1*Lamda|1 - H2*L|
|00001140| 61 6d 64 61 32 0d 07 62 | 1c 48 32 3d 48 32 2a 4c |amda2..b|.H2=H2*L|
|00001150| 61 6d 64 61 31 20 2b 20 | 48 31 2a 4c 61 6d 64 61 |amda1 + |H1*Lamda|
|00001160| 32 0d 07 6c 08 48 31 3d | 54 0d 07 76 0c 5a 31 3d |2..l.H1=|T..v.Z1=|
|00001170| 5a 31 2b 48 31 0d 07 80 | 0c 5a 32 3d 5a 32 2b 48 |Z1+H1...|.Z2=Z2+H|
|00001180| 32 0d 07 8a 05 20 0d 07 | 94 05 f5 0d 07 9e 05 20 |2.... ..|....... |
|00001190| 0d 07 a8 19 64 69 76 65 | 72 3d 30 20 3a 20 49 74 |....dive|r=0 : It|
|000011a0| 65 72 3d 49 74 65 72 2b | 31 0d 07 b2 09 f2 65 76 |er=Iter+|1.....ev|
|000011b0| 61 6c 0d 07 bc 0b f2 6f | 75 74 70 75 74 0d 07 c6 |al.....o|utput...|
|000011c0| 0f e7 20 61 62 6f 72 74 | 20 8c 20 e1 0d 07 d0 05 |.. abort| . .....|
|000011d0| 20 0d 07 da 18 f4 20 43 | 6f 6e 76 65 72 67 65 6e | ..... C|onvergen|
|000011e0| 63 65 20 43 68 65 63 6b | 73 0d 07 e4 04 0d 07 ee |ce Check|s.......|
|000011f0| 33 48 6e 6f 72 6d 3d 94 | 28 48 31 29 2b 94 28 48 |3Hnorm=.|(H1)+.(H|
|00001200| 32 29 20 3a 20 f4 20 73 | 68 69 66 74 20 69 6e 20 |2) : . s|hift in |
|00001210| 69 74 65 72 61 74 65 20 | 70 6f 73 69 74 69 6f 6e |iterate |position|
|00001220| 0d 07 f8 32 5a 6e 6f 72 | 6d 3d 94 28 5a 31 29 2b |...2Znor|m=.(Z1)+|
|00001230| 94 28 5a 32 29 20 3a 20 | f4 20 63 75 72 72 65 6e |.(Z2) : |. curren|
|00001240| 74 20 69 74 65 72 61 74 | 65 20 70 6f 73 69 74 69 |t iterat|e positi|
|00001250| 6f 6e 0d 08 02 34 46 6e | 6f 72 6d 3d 94 28 46 31 |on...4Fn|orm=.(F1|
|00001260| 29 2b 94 28 46 32 29 20 | 3a 20 f4 20 20 20 20 64 |)+.(F2) |: . d|
|00001270| 65 66 6c 61 74 65 64 20 | 66 75 6e 63 74 69 6f 6e |eflated |function|
|00001280| 20 76 61 6c 75 65 0d 08 | 0c 34 55 6e 6f 72 6d 3d | value..|.4Unorm=|
|00001290| 94 28 55 31 29 2b 94 28 | 55 32 29 20 3a 20 f4 20 |.(U1)+.(|U2) : . |
|000012a0| 75 6e 2d 64 65 66 6c 61 | 74 65 64 20 66 75 6e 63 |un-defla|ted func|
|000012b0| 74 69 6f 6e 20 76 61 6c | 75 65 0d 08 16 37 47 6e |tion val|ue...7Gn|
|000012c0| 6f 72 6d 3d 94 28 4f 6c | 64 66 31 29 2b 94 28 4f |orm=.(Ol|df1)+.(O|
|000012d0| 6c 64 66 32 29 20 3a 20 | f4 20 70 72 65 76 69 6f |ldf2) : |. previo|
|000012e0| 75 73 20 66 75 6e 63 74 | 69 6f 6e 20 76 61 6c 75 |us funct|ion valu|
|000012f0| 65 0d 08 20 31 63 6f 6e | 76 3d 28 48 6e 6f 72 6d |e.. 1con|v=(Hnorm|
|00001300| 3c 45 31 2a 5a 6e 6f 72 | 6d 29 20 84 20 28 46 6e |<E1*Znor|m) . (Fn|
|00001310| 6f 72 6d 3c 45 32 20 80 | 20 55 6e 6f 72 6d 3c 45 |orm<E2 .| Unorm<E|
|00001320| 32 29 0d 08 2a 04 0d 08 | 34 1e e7 20 46 6e 6f 72 |2)..*...|4.. Fnor|
|00001330| 6d 3e 31 30 2a 47 6e 6f | 72 6d 20 84 20 65 72 72 |m>10*Gno|rm . err|
|00001340| 3d 32 20 8c 0d 08 3e 0c | e7 20 63 6f 6e 76 20 8c |=2 ...>.|. conv .|
|00001350| 0d 08 48 0c 5a 31 3d 5a | 31 2d 48 31 0d 08 52 0c |..H.Z1=Z|1-H1..R.|
|00001360| 5a 32 3d 5a 32 2d 48 32 | 0d 08 5c 05 cc 0d 08 66 |Z2=Z2-H2|..\....f|
|00001370| 0b 64 69 76 65 72 3d 31 | 0d 08 70 0b 48 31 3d 48 |.diver=1|..p.H1=H|
|00001380| 31 2f 32 0d 08 7a 0b 48 | 32 3d 48 32 2f 32 0d 08 |1/2..z.H|2=H2/2..|
|00001390| 84 45 f4 20 6e 65 78 74 | 20 6c 69 6e 65 20 66 72 |.E. next| line fr|
|000013a0| 6f 6d 20 43 6f 6e 74 65 | 20 26 20 64 65 20 42 6f |om Conte| & de Bo|
|000013b0| 6f 72 20 62 75 74 20 70 | 72 6f 62 6c 65 6d 73 20 |or but p|roblems |
|000013c0| 6f 6e 20 68 69 67 68 2d | 6f 72 64 65 72 20 70 6f |on high-|order po|
|000013d0| 6c 79 73 0d 08 8e 25 f4 | 20 4c 61 6d 64 61 31 3d |lys...%.| Lamda1=|
|000013e0| 4c 61 6d 64 61 31 2f 32 | 3a 4c 61 6d 64 61 32 3d |Lamda1/2|:Lamda2=|
|000013f0| 4c 61 6d 64 61 32 2f 32 | 0d 08 98 0c 5a 31 3d 5a |Lamda2/2|....Z1=Z|
|00001400| 31 2d 48 31 0d 08 a2 0c | 5a 32 3d 5a 32 2d 48 32 |1-H1....|Z2=Z2-H2|
|00001410| 0d 08 ac 0a 63 66 6c 67 | 3d 30 0d 08 b6 09 65 72 |....cflg|=0....er|
|00001420| 72 3d 30 0d 08 c0 05 cd | 0d 08 ca 05 cd 0d 08 d4 |r=0.....|........|
|00001430| 04 0d 08 de 13 fd 20 64 | 69 76 65 72 3d 30 20 84 |...... d|iver=0 .|
|00001440| 20 65 72 72 0d 08 e8 05 | 20 0d 08 f2 19 fd 20 28 | err....| ..... (|
|00001450| 63 6f 6e 76 20 80 20 63 | 66 6c 67 29 20 84 20 65 |conv . c|flg) . e|
|00001460| 72 72 0d 08 fc 05 20 0d | 09 06 0b fd 20 65 72 72 |rr.... .|.... err|
|00001470| 3d 30 0d 09 10 05 20 0d | 09 1a 10 5a 72 31 28 72 |=0.... .|...Zr1(r|
|00001480| 6f 6f 74 29 3d 5a 31 0d | 09 24 10 5a 72 32 28 72 |oot)=Z1.|.$.Zr2(r|
|00001490| 6f 6f 74 29 3d 5a 32 0d | 09 2e 05 20 0d 09 38 1a |oot)=Z2.|... ..8.|
|000014a0| fb 20 32 20 3a 20 f1 20 | 27 22 20 20 43 6f 6e 76 |. 2 : . |'" Conv|
|000014b0| 65 72 67 65 64 22 0d 09 | 42 07 fb 20 30 0d 09 4c |erged"..|B.. 0..L|
|000014c0| 15 f1 20 27 22 20 50 72 | 65 73 73 20 61 20 6b 65 |.. '" Pr|ess a ke|
|000014d0| 79 22 0d 09 56 0f e7 20 | 72 6f 6f 74 3c 4e 72 20 |y"..V.. |root<Nr |
|000014e0| 8c 0d 09 60 15 f1 20 22 | 20 74 6f 20 63 6f 6e 74 |...`.. "| to cont|
|000014f0| 69 6e 75 65 22 3b 0d 09 | 6a 05 cc 0d 09 74 15 f1 |inue";..|j....t..|
|00001500| 20 22 20 74 6f 20 45 6e | 64 20 50 72 6f 67 22 3b | " to En|d Prog";|
|00001510| 0d 09 7e 0d 40 25 3d 31 | 30 3a ef 32 30 0d 09 88 |..~.@%=1|0:.20...|
|00001520| 05 cd 0d 09 92 0b 2a 46 | 58 31 35 2c 31 0d 09 9c |......*F|X15,1...|
|00001530| 11 d4 20 31 2c 2d 31 35 | 2c 31 35 30 2c 35 0d 09 |.. 1,-15|,150,5..|
|00001540| a6 0b 47 3d a5 20 3a 20 | db 0d 09 b0 05 20 0d 09 |..G=. : |..... ..|
|00001550| ba 0a ed 20 72 6f 6f 74 | 0d 09 c4 05 20 0d 09 ce |... root|.... ...|
|00001560| 05 e1 0d 09 d8 04 0d 09 | e2 04 0d 09 ec 0b dd 20 |........|....... |
|00001570| f2 65 76 61 6c 0d 09 f6 | 1a e7 20 46 24 3d 22 46 |.eval...|.. F$="F|
|00001580| 22 20 f2 66 75 6e 63 20 | 8b 20 f2 70 6f 6c 79 0d |" .func |. .poly.|
|00001590| 0a 00 0c f2 73 75 70 72 | 65 73 73 0d 0a 0a 14 e7 |....supr|ess.....|
|000015a0| 20 64 65 62 75 67 20 8c | 20 f2 64 65 62 75 67 0d | debug .| .debug.|
|000015b0| 0a 14 05 e1 0d 0a 1e 05 | 20 0d 0a 28 04 0d 0a 32 |........| ..(...2|
|000015c0| 0c dd 20 f2 64 65 62 75 | 67 0d 0a 3c 0e c8 8e 20 |.. .debu|g..<... |
|000015d0| 64 65 62 75 67 20 ca 0d | 0a 46 09 c9 20 31 20 3a |debug ..|.F.. 1 :|
|000015e0| 0d 0a 50 06 ef 32 0d 0a | 5a 5f f1 20 22 5a 31 3d |..P..2..|Z_. "Z1=|
|000015f0| 22 3b 5a 31 27 22 5a 32 | 3d 22 3b 5a 32 27 22 55 |";Z1'"Z2|=";Z2'"U|
|00001600| 31 3d 22 3b 55 31 27 22 | 55 32 3d 22 3b 55 32 27 |1=";U1'"|U2=";U2'|
|00001610| 22 48 31 3d 22 3b 48 31 | 27 22 48 32 3d 22 3b 48 |"H1=";H1|'"H2=";H|
|00001620| 32 27 22 4c 61 6d 64 61 | 31 3d 22 3b 4c 61 6d 64 |2'"Lamda|1=";Lamd|
|00001630| 61 31 27 22 4c 61 6d 64 | 61 32 3d 22 3b 4c 61 6d |a1'"Lamd|a2=";Lam|
|00001640| 64 61 32 27 27 0d 0a 64 | 06 ef 33 0d 0a 6e 13 ff |da2''..d|..3..n..|
|00001650| 28 22 46 58 31 35 2c 31 | 22 29 3a 47 3d a5 0d 0a |("FX15,1|"):G=...|
|00001660| 78 09 c9 20 32 20 3a 0d | 0a 82 25 e3 20 72 6f 77 |x.. 2 :.|..%. row|
|00001670| 3d 32 36 20 b8 20 33 30 | 3a f1 20 8a 30 2c 72 6f |=26 . 30|:. .0,ro|
|00001680| 77 29 3b 42 6c 61 6e 6b | 24 3b 3a ed 0d 0a 8c 38 |w);Blank|$;:....8|
|00001690| f1 20 8a 30 2c 32 36 29 | 3b 22 5a 31 3d 22 3b 5a |. .0,26)|;"Z1=";Z|
|000016a0| 31 27 22 5a 32 3d 22 3b | 5a 32 27 22 55 31 3d 22 |1'"Z2=";|Z2'"U1="|
|000016b0| 3b 55 31 27 22 55 32 3d | 22 3b 55 32 3a 2a 46 58 |;U1'"U2=|";U2:*FX|
|000016c0| 31 35 2c 31 0d 0a 96 07 | 47 3d a5 0d 0a a0 25 e3 |15,1....|G=....%.|
|000016d0| 20 72 6f 77 3d 32 36 20 | b8 20 33 30 3a f1 20 8a | row=26 |. 30:. .|
|000016e0| 30 2c 72 6f 77 29 3b 42 | 6c 61 6e 6b 24 3b 3a ed |0,row);B|lank$;:.|
|000016f0| 0d 0a aa 0c f1 20 8a 30 | 2c 30 29 3b 0d 0a b4 05 |..... .0|,0);....|
|00001700| cb 0d 0a be 05 e1 0d 0a | c8 04 0d 0a d2 04 0d 0a |........|........|
|00001710| dc 35 f4 20 45 76 61 6c | 75 61 74 65 20 61 20 72 |.5. Eval|uate a r|
|00001720| 65 61 6c 20 70 6f 6c 79 | 6e 6f 6d 69 61 6c 20 66 |eal poly|nomial f|
|00001730| 6f 72 20 63 6f 6d 70 6c | 65 78 20 61 72 67 75 6d |or compl|ex argum|
|00001740| 65 6e 74 0d 0a e6 33 f4 | 20 75 73 69 6e 67 20 74 |ent...3.| using t|
|00001750| 68 65 20 70 6f 6c 61 72 | 20 65 78 70 61 6e 73 69 |he polar| expansi|
|00001760| 6f 6e 20 61 6c 67 6f 72 | 69 74 68 6d 20 77 68 69 |on algor|ithm whi|
|00001770| 63 68 20 68 61 73 0d 0a | f0 2b f4 20 67 6f 6f 64 |ch has..|.+. good|
|00001780| 20 61 63 63 75 72 61 63 | 79 20 66 6f 72 20 7c 5a | accurac|y for |Z|
|00001790| 7c 20 63 6c 6f 73 65 20 | 74 6f 20 75 6e 69 74 79 || close |to unity|
|000017a0| 2e 0d 0a fa 04 0d 0b 04 | 0b dd 20 f2 70 6f 6c 79 |........|.. .poly|
|000017b0| 0d 0b 0e 05 20 0d 0b 18 | 09 65 72 72 3d 30 0d 0b |.... ...|.err=0..|
|000017c0| 22 19 55 31 3d 43 28 4e | 72 29 20 3a 20 55 32 3d |".U1=C(N|r) : U2=|
|000017d0| 43 28 4e 72 2d 31 29 0d | 0b 2c 20 54 77 6f 78 3d |C(Nr-1).|., Twox=|
|000017e0| 32 2a 5a 31 20 3a 20 52 | 61 64 32 3d 5a 31 2a 5a |2*Z1 : R|ad2=Z1*Z|
|000017f0| 31 2b 5a 32 2a 5a 32 0d | 0b 36 1b e7 20 52 61 64 |1+Z2*Z2.|.6.. Rad|
|00001800| 32 3e 4d 61 78 6d 6f 64 | 20 8c 20 65 72 72 3d 32 |2>Maxmod| . err=2|
|00001810| 3a e1 0d 0b 40 04 0d 0b | 4a 15 e3 20 4a 3d 4e 72 |:...@...|J.. J=Nr|
|00001820| 2d 32 20 b8 20 30 20 88 | 20 2d 31 0d 0b 54 13 54 |-2 . 0 .| -1..T.T|
|00001830| 31 3d 55 32 20 2b 20 54 | 77 6f 78 2a 55 31 0d 0b |1=U2 + T|wox*U1..|
|00001840| 5e 13 55 32 3d 43 28 4a | 29 2d 52 61 64 32 2a 55 |^.U2=C(J|)-Rad2*U|
|00001850| 31 0d 0b 68 09 55 31 3d | 54 31 0d 0b 72 07 ed 20 |1..h.U1=|T1..r.. |
|00001860| 4a 0d 0b 7c 0f 54 31 3d | 5a 31 2a 55 31 2b 55 32 |J..|.T1=|Z1*U1+U2|
|00001870| 0d 0b 86 0c 55 32 3d 5a | 32 2a 55 31 0d 0b 90 09 |....U2=Z|2*U1....|
|00001880| 55 31 3d 54 31 0d 0b 9a | 04 0d 0b a4 1d e7 20 94 |U1=T1...|...... .|
|00001890| 28 55 31 29 2b 94 28 55 | 32 29 3e 45 31 39 20 8c |(U1)+.(U|2)>E19 .|
|000018a0| 20 65 72 72 3d 32 0d 0b | ae 04 0d 0b b8 05 e1 0d | err=2..|........|
|000018b0| 0b c2 05 20 0d 0b cc 04 | 0d 0b d6 46 f4 20 53 75 |... ....|...F. Su|
|000018c0| 70 70 72 65 73 73 20 70 | 72 65 76 69 6f 75 73 6c |ppress p|reviousl|
|000018d0| 79 20 64 69 73 63 6f 76 | 65 72 65 64 20 72 6f 6f |y discov|ered roo|
|000018e0| 74 73 20 74 6f 20 67 69 | 76 65 20 64 65 66 6c 61 |ts to gi|ve defla|
|000018f0| 74 65 64 20 66 75 6e 63 | 74 69 6f 6e 20 46 0d 0b |ted func|tion F..|
|00001900| e0 04 0d 0b ea 0e dd 20 | f2 73 75 70 72 65 73 73 |....... |.supress|
|00001910| 0d 0b f4 05 20 0d 0b fe | 0f 46 31 3d 55 31 3a 50 |.... ...|.F1=U1:P|
|00001920| 72 31 3d 31 0d 0c 08 0f | 46 32 3d 55 32 3a 50 72 |r1=1....|F2=U2:Pr|
|00001930| 32 3d 30 0d 0c 12 05 20 | 0d 0c 1c 14 e7 20 65 72 |2=0.... |..... er|
|00001940| 72 20 84 20 72 6f 6f 74 | 3d 31 20 e1 0d 0c 26 05 |r . root|=1 ...&.|
|00001950| 20 0d 0c 30 12 e3 20 4a | 3d 31 20 b8 20 72 6f 6f | ..0.. J|=1 . roo|
|00001960| 74 2d 31 0d 0c 3a 10 44 | 31 3d 5a 31 2d 5a 72 31 |t-1..:.D|1=Z1-Zr1|
|00001970| 28 4a 29 0d 0c 44 10 44 | 32 3d 5a 32 2d 5a 72 32 |(J)..D.D|2=Z2-Zr2|
|00001980| 28 4a 29 0d 0c 4e 15 54 | 72 31 3d 50 72 31 2a 44 |(J)..N.T|r1=Pr1*D|
|00001990| 31 2d 50 72 32 2a 44 32 | 0d 0c 58 15 50 72 32 3d |1-Pr2*D2|..X.Pr2=|
|000019a0| 50 72 32 2a 44 31 2b 50 | 72 31 2a 44 32 0d 0c 62 |Pr2*D1+P|r1*D2..b|
|000019b0| 0b 50 72 31 3d 54 72 31 | 0d 0c 6c 07 ed 20 4a 0d |.Pr1=Tr1|..l.. J.|
|000019c0| 0c 76 19 50 72 6d 3d 28 | 50 72 31 2a 50 72 31 2b |.v.Prm=(|Pr1*Pr1+|
|000019d0| 50 72 32 2a 50 72 32 29 | 0d 0c 80 0d e7 20 50 72 |Pr2*Pr2)|..... Pr|
|000019e0| 6d 3d 30 20 8c 0d 0c 8a | 09 65 72 72 3d 32 0d 0c |m=0 ....|.err=2..|
|000019f0| 94 05 cc 0d 0c 9e 14 e7 | 20 50 72 6d 3c 45 33 20 |........| Prm<E3 |
|00001a00| 8c 20 65 72 72 3d 31 0d | 0c a8 05 cd 0d 0c b2 0d |. err=1.|........|
|00001a10| e7 20 65 72 72 3d 30 20 | 8c 0d 0c bc 1a 46 31 3d |. err=0 |.....F1=|
|00001a20| 28 55 31 2a 50 72 31 2b | 55 32 2a 50 72 32 29 2f |(U1*Pr1+|U2*Pr2)/|
|00001a30| 50 72 6d 0d 0c c6 1a 46 | 32 3d 28 55 32 2a 50 72 |Prm....F|2=(U2*Pr|
|00001a40| 31 2d 55 31 2a 50 72 32 | 29 2f 50 72 6d 0d 0c d0 |1-U1*Pr2|)/Prm...|
|00001a50| 05 cd 0d 0c da 04 0d 0c | e4 05 e1 0d 0c ee 04 0d |........|........|
|00001a60| 0c f8 04 0d 0d 02 41 f4 | 20 53 75 62 72 6f 75 74 |......A.| Subrout|
|00001a70| 69 6e 65 20 66 6f 72 20 | 63 6f 6d 70 6c 65 78 20 |ine for |complex |
|00001a80| 73 71 75 61 72 65 2d 72 | 6f 6f 74 20 28 69 6e 63 |square-r|oot (inc|
|00001a90| 6c 75 64 65 73 20 76 61 | 72 69 6f 75 73 20 74 72 |ludes va|rious tr|
|00001aa0| 61 70 73 29 0d 0d 0c 04 | 0d 0d 16 0b dd 20 f2 73 |aps)....|..... .s|
|00001ab0| 71 72 74 0d 0d 20 05 20 | 0d 0d 2a 11 53 71 31 3d |qrt.. . |..*.Sq1=|
|00001ac0| 30 20 3a 20 53 71 32 3d | 30 0d 0d 34 19 41 62 31 |0 : Sq2=|0..4.Ab1|
|00001ad0| 3d 94 28 42 31 29 20 3a | 20 41 62 32 3d 94 28 42 |=.(B1) :| Ab2=.(B|
|00001ae0| 32 29 0d 0d 3e 05 20 0d | 0d 48 14 e7 20 41 62 31 |2)..>. .|.H.. Ab1|
|00001af0| 2b 41 62 32 20 3c 20 45 | 34 20 e1 0d 0d 52 1e e7 |+Ab2 < E|4 ...R..|
|00001b00| 20 42 32 3d 30 20 80 20 | 42 31 3e 3d 30 20 53 71 | B2=0 . |B1>=0 Sq|
|00001b10| 31 3d b6 28 42 31 29 3a | e1 0d 0d 5c 1d e7 20 42 |1=.(B1):|...\.. B|
|00001b20| 32 3d 30 20 53 71 32 3d | b6 28 2d 42 31 2a 63 6d |2=0 Sq2=|.(-B1*cm|
|00001b30| 70 6c 78 29 3a e1 0d 0d | 66 05 20 0d 0d 70 13 e7 |plx):...|f. ..p..|
|00001b40| 20 41 62 31 2b 41 62 32 | 3e 45 31 39 20 8c 0d 0d | Ab1+Ab2|>E19 ...|
|00001b50| 7a 2d e7 20 41 62 31 3e | 41 62 32 20 8c 20 56 76 |z-. Ab1>|Ab2 . Vv|
|00001b60| 3d 41 62 31 3a 57 77 3d | 41 62 32 20 8b 20 56 76 |=Ab1:Ww=|Ab2 . Vv|
|00001b70| 3d 41 62 32 3a 57 77 3d | 41 62 31 0d 0d 84 0c 58 |=Ab2:Ww=|Ab1....X|
|00001b80| 78 3d 57 77 2f 56 76 0d | 0d 8e 15 53 62 30 3d 56 |x=Ww/Vv.|...Sb0=V|
|00001b90| 76 2a b6 28 31 2b 58 78 | 2a 58 78 29 0d 0d 98 05 |v*.(1+Xx|*Xx)....|
|00001ba0| cc 0d 0d a2 18 53 62 30 | 3d b6 28 42 31 2a 42 31 |.....Sb0|=.(B1*B1|
|00001bb0| 20 2b 20 42 32 2a 42 32 | 29 0d 0d ac 05 cd 0d 0d | + B2*B2|).......|
|00001bc0| b6 21 53 71 31 3d b4 28 | 42 32 29 20 2a 20 b6 28 |.!Sq1=.(|B2) * .(|
|00001bd0| 20 2e 35 2a 94 28 42 31 | 2b 53 62 30 29 20 29 0d | .5*.(B1|+Sb0) ).|
|00001be0| 0d c0 0d e7 20 53 71 31 | 3d 30 20 8c 0d 0d ca 10 |.... Sq1|=0 .....|
|00001bf0| 53 71 32 3d b6 28 94 28 | 42 31 29 29 0d 0d d4 05 |Sq2=.(.(|B1))....|
|00001c00| cc 0d 0d de 14 53 71 32 | 3d 94 28 2e 35 2a 42 32 |.....Sq2|=.(.5*B2|
|00001c10| 2f 53 71 31 29 0d 0d e8 | 05 cd 0d 0d f2 05 20 0d |/Sq1)...|...... .|
|00001c20| 0d fc 05 e1 0d 0e 06 05 | 20 0d 0e 10 04 0d 0e 1a |........| .......|
|00001c30| 0c dd 20 f2 63 6f 65 66 | 66 0d 0e 24 05 20 0d 0e |.. .coef|f..$. ..|
|00001c40| 2e 13 e3 20 4a 3d 4e 72 | 20 b8 20 30 20 88 20 2d |... J=Nr| . 0 . -|
|00001c50| 31 0d 0e 38 0a f3 20 43 | 28 4a 29 0d 0e 42 07 ed |1..8.. C|(J)..B..|
|00001c60| 20 4a 0d 0e 4c 05 20 0d | 0e 56 05 e1 0d 0e 60 05 | J..L. .|.V....`.|
|00001c70| 20 0d 0e 6a 04 0d 0e 74 | 0d dd 20 f2 6f 75 74 70 | ..j...t|.. .outp|
|00001c80| 75 74 0d 0e 7e 05 20 0d | 0e 88 15 ef 33 30 20 3a |ut..~. .|....30 :|
|00001c90| 20 fb 20 32 20 3a 20 40 | 25 3d 31 30 0d 0e 92 17 | . 2 : @|%=10....|
|00001ca0| f1 20 27 22 20 20 20 20 | 52 4f 4f 54 20 22 3b 72 |. '" |ROOT ";r|
|00001cb0| 6f 6f 74 0d 0e 9c 1c f1 | 20 27 3b 22 20 20 49 74 |oot.....| ';" It|
|00001cc0| 65 72 61 74 65 20 22 3b | 49 74 65 72 27 27 27 0d |erate ";|Iter'''.|
|00001cd0| 0e a6 18 4d 6f 64 3d b6 | 28 5a 31 2a 5a 31 20 2b |...Mod=.|(Z1*Z1 +|
|00001ce0| 20 5a 32 2a 5a 32 29 0d | 0e b0 13 e7 20 5a 31 3d | Z2*Z2).|.... Z1=|
|00001cf0| 30 20 84 20 5a 32 3d 30 | 20 8c 0d 0e ba 19 e7 20 |0 . Z2=0| ...... |
|00001d00| 5a 31 3d 30 20 54 3d 31 | 38 30 2d 39 30 2a b4 28 |Z1=0 T=1|80-90*.(|
|00001d10| 5a 32 29 0d 0e c4 19 e7 | 20 5a 32 3d 30 20 54 3d |Z2).....| Z2=0 T=|
|00001d20| 39 30 20 2d 39 30 2a b4 | 28 5a 31 29 0d 0e ce 05 |90 -90*.|(Z1)....|
|00001d30| cc 0d 0e d8 13 54 3d 9d | 28 20 99 28 5a 32 2f 5a |.....T=.|( .(Z2/Z|
|00001d40| 31 29 20 29 0d 0e e2 0c | e7 20 5a 31 3c 30 20 8c |1) )....|. Z1<0 .|
|00001d50| 0d 0e ec 0b 54 3d 31 38 | 30 2b 54 0d 0e f6 05 cc |....T=18|0+T.....|
|00001d60| 0d 0f 00 1d 54 3d 54 2b | 33 36 30 20 3a 20 e7 20 |....T=T+|360 : . |
|00001d70| 54 3e 33 36 30 20 54 3d | 54 2d 33 36 30 0d 0f 0a |T>360 T=|T-360...|
|00001d80| 05 cd 0d 0f 14 05 cd 0d | 0f 1e 0d 40 25 3d 26 32 |........|...@%=&2|
|00001d90| 30 38 30 44 0d 0f 28 17 | f1 20 22 20 20 20 20 52 |080D..(.|. " R|
|00001da0| 65 28 5a 29 22 20 3a 20 | fb 20 33 0d 0f 32 12 e7 |e(Z)" : |. 3..2..|
|00001db0| 20 5a 31 3e 3d 30 20 f1 | 20 22 20 22 3b 0d 0f 3c | Z1>=0 .| " ";..<|
|00001dc0| 12 f1 20 22 20 22 3b 5a | 31 20 3a 20 fb 20 32 0d |.. " ";Z|1 : . 2.|
|00001dd0| 0f 46 17 f1 20 22 20 20 | 20 20 49 6d 28 5a 29 22 |.F.. " | Im(Z)"|
|00001de0| 20 3a 20 fb 20 33 0d 0f | 50 12 e7 20 5a 32 3e 3d | : . 3..|P.. Z2>=|
|00001df0| 30 20 f1 20 22 20 22 3b | 0d 0f 5a 15 f1 20 22 20 |0 . " ";|..Z.. " |
|00001e00| 22 3b 5a 32 20 27 27 20 | 3a 20 fb 20 32 0d 0f 64 |";Z2 '' |: . 2..d|
|00001e10| 1a f1 20 22 20 20 4d 6f | 64 75 6c 75 73 28 5a 29 |.. " Mo|dulus(Z)|
|00001e20| 22 20 3a 20 fb 20 33 0d | 0f 6e 14 f1 20 22 20 20 |" : . 3.|.n.. " |
|00001e30| 22 3b 4d 6f 64 20 3a 20 | fb 20 32 0d 0f 78 1a f1 |";Mod : |. 2..x..|
|00001e40| 20 22 20 20 41 72 67 28 | 5a 29 20 64 65 67 22 20 | " Arg(|Z) deg" |
|00001e50| 3a 20 fb 20 33 0d 0f 82 | 1a f1 20 22 20 20 20 20 |: . 3...|.. " |
|00001e60| 20 20 20 20 20 20 20 20 | 20 22 3a 20 ef 31 31 0d | | ": .11.|
|00001e70| 0f 8c 12 f1 20 22 20 20 | 22 3b 54 20 3a 20 fb 20 |.... " |";T : . |
|00001e80| 32 0d 0f 96 08 f1 20 27 | 27 0d 0f a0 0d 40 25 3d |2..... '|'....@%=|
|00001e90| 26 31 30 36 30 44 0d 0f | aa 17 f1 20 22 20 20 20 |&1060D..|... " |
|00001ea0| 20 52 65 28 46 29 22 20 | 3a 20 fb 20 33 0d 0f b4 | Re(F)" |: . 3...|
|00001eb0| 12 e7 20 55 31 3e 3d 30 | 20 f1 20 22 20 22 3b 0d |.. U1>=0| . " ";.|
|00001ec0| 0f be 12 f1 20 22 20 22 | 3b 55 31 20 3a 20 fb 20 |.... " "|;U1 : . |
|00001ed0| 32 0d 0f c8 17 f1 20 22 | 20 20 20 20 49 6d 28 46 |2..... "| Im(F|
|00001ee0| 29 22 20 3a 20 fb 20 33 | 0d 0f d2 12 e7 20 55 32 |)" : . 3|..... U2|
|00001ef0| 3e 3d 30 20 f1 20 22 20 | 22 3b 0d 0f dc 10 f1 20 |>=0 . " |";..... |
|00001f00| 22 20 22 3b 55 32 20 3a | 20 f1 0d 0f e6 05 20 0d |" ";U2 :| ..... .|
|00001f10| 0f f0 1a e7 20 49 74 65 | 72 3d 31 20 58 70 3d 45 |.... Ite|r=1 Xp=E|
|00001f20| 36 20 3a 20 59 70 3d 45 | 36 0d 0f fa 20 c8 92 20 |6 : Yp=E|6... .. |
|00001f30| 58 73 2a 28 58 70 2d 58 | 6d 69 6e 29 2c 59 73 2a |Xs*(Xp-X|min),Ys*|
|00001f40| 28 59 70 2d 59 6d 69 6e | 29 0d 10 04 1c e7 20 94 |(Yp-Ymin|)..... .|
|00001f50| 28 58 73 2a 5a 31 29 2b | 94 28 59 73 2a 5a 32 29 |(Xs*Z1)+|.(Ys*Z2)|
|00001f60| 3e 45 36 20 8c 0d 10 0e | 2c 58 70 3d 45 36 3a 59 |>E6 ....|,Xp=E6:Y|
|00001f70| 70 3d 45 36 3a c8 92 20 | 58 73 2a 28 58 70 2d 58 |p=E6:.. |Xs*(Xp-X|
|00001f80| 6d 69 6e 29 2c 59 73 2a | 28 59 70 2d 59 6d 69 6e |min),Ys*|(Yp-Ymin|
|00001f90| 29 0d 10 18 05 cc 0d 10 | 22 20 c8 92 20 58 73 2a |).......|" .. Xs*|
|00001fa0| 28 5a 31 2d 58 6d 69 6e | 29 2c 59 73 2a 28 5a 32 |(Z1-Xmin|),Ys*(Z2|
|00001fb0| 2d 59 6d 69 6e 29 0d 10 | 2c 09 58 70 3d 5a 31 0d |-Ymin)..|,.Xp=Z1.|
|00001fc0| 10 36 09 59 70 3d 5a 32 | 0d 10 40 05 cd 0d 10 4a |.6.Yp=Z2|..@....J|
|00001fd0| 05 20 0d 10 54 18 e7 20 | 49 74 65 72 20 83 20 4d |. ..T.. |Iter . M|
|00001fe0| 61 78 69 74 20 3d 20 30 | 20 8c 0d 10 5e 07 fb 20 |axit = 0| ...^.. |
|00001ff0| 32 0d 10 68 22 f1 27 22 | 20 49 74 65 72 20 63 6f |2..h".'"| Iter co|
|00002000| 75 6e 74 22 27 22 20 75 | 70 20 74 6f 20 4d 61 78 |unt"'" u|p to Max|
|00002010| 69 74 22 0d 10 72 07 fb | 20 30 0d 10 7c 16 f5 20 |it"..r..| 0..|.. |
|00002020| 3a 20 f1 20 8a 30 2c 33 | 30 29 3b 20 3a 20 ef 37 |: . .0,3|0); : .7|
|00002030| 0d 10 86 0b 2a 46 58 31 | 35 2c 31 0d 10 90 2c f1 |....*FX1|5,1...,.|
|00002040| 22 20 45 6e 74 65 72 20 | 53 28 74 6f 70 29 22 27 |" Enter |S(top)"'|
|00002050| 22 20 20 20 20 6f 72 20 | 43 28 6f 6e 74 29 22 3b |" or |C(ont)";|
|00002060| 20 3a 20 51 24 3d be 0d | 10 9a 27 fd 20 51 24 3d | : Q$=..|..'. Q$=|
|00002070| 22 53 22 20 84 20 51 24 | 3d 22 43 22 20 84 20 51 |"S" . Q$|="C" . Q|
|00002080| 24 3d 22 73 22 20 84 20 | 51 24 3d 22 63 22 0d 10 |$="s" . |Q$="c"..|
|00002090| a4 17 e7 20 51 24 3d 22 | 53 22 20 84 20 51 24 3d |... Q$="|S" . Q$=|
|000020a0| 22 73 22 20 8c 0d 10 ae | 0b 61 62 6f 72 74 3d 31 |"s" ....|.abort=1|
|000020b0| 0d 10 b8 05 cc 0d 10 c2 | 25 e3 20 72 6f 77 3d 32 |........|%. row=2|
|000020c0| 36 20 b8 20 33 30 3a f1 | 20 8a 30 2c 72 6f 77 29 |6 . 30:.| .0,row)|
|000020d0| 3b 42 6c 61 6e 6b 24 3b | 3a ed 0d 10 cc 05 cd 0d |;Blank$;|:.......|
|000020e0| 10 d6 05 cd 0d 10 e0 04 | 0d 10 ea 05 e1 0d 10 f4 |........|........|
|000020f0| 05 20 0d 10 fe 04 0d 11 | 08 12 dd 20 f2 73 63 72 |. ......|... .scr|
|00002100| 65 65 6e 5f 69 6e 69 74 | 0d 11 12 05 20 0d 11 1c |een_init|.... ...|
|00002110| 15 53 63 6d 3d 38 30 30 | 20 3a 20 61 62 6f 72 74 |.Scm=800| : abort|
|00002120| 3d 30 0d 11 26 16 58 73 | 3d 53 63 6d 2f 28 58 6d |=0..&.Xs|=Scm/(Xm|
|00002130| 61 78 2d 58 6d 69 6e 29 | 0d 11 30 16 59 73 3d 53 |ax-Xmin)|..0.Ys=S|
|00002140| 63 6d 2f 28 59 6d 61 78 | 2d 59 6d 69 6e 29 0d 11 |cm/(Ymax|-Ymin)..|
|00002150| 3a 11 da 20 3a 20 fb 20 | 31 32 39 20 3a 20 db 0d |:.. : . |129 : ..|
|00002160| 11 44 13 ec 20 30 2c 30 | 20 3a 20 df 20 53 63 6d |.D.. 0,0| : . Scm|
|00002170| 2c 30 0d 11 4e 0e df 20 | 53 63 6d 2c 31 30 32 33 |,0..N.. |Scm,1023|
|00002180| 0d 11 58 14 df 20 30 2c | 31 30 32 33 20 3a 20 df |..X.. 0,|1023 : .|
|00002190| 20 30 2c 30 0d 11 62 17 | ec 20 30 2c 39 35 30 20 | 0,0..b.|. 0,950 |
|000021a0| 3a 20 df 20 53 63 6d 2c | 39 35 30 0d 11 6c 17 ec |: . Scm,|950..l..|
|000021b0| 20 30 2c 53 63 6d 20 3a | 20 df 20 53 63 6d 2c 53 | 0,Scm :| . Scm,S|
|000021c0| 63 6d 0d 11 76 05 20 0d | 11 80 1b ef 20 35 20 3a |cm..v. .|.... 5 :|
|000021d0| 20 ec 20 36 30 2c 31 30 | 30 30 20 3a 20 e6 20 30 | . 60,10|00 : . 0|
|000021e0| 2c 32 0d 11 8a 1e f1 20 | 22 50 72 6f 67 72 61 6d |,2..... |"Program|
|000021f0| 20 5a 65 72 46 75 6e 63 | 20 3a 20 22 3b 66 6e 24 | ZerFunc| : ";fn$|
|00002200| 0d 11 94 2f ec 20 38 30 | 2c 38 36 35 20 3a 20 f1 |.../. 80|,865 : .|
|00002210| 22 58 6d 69 6e 3d 22 3b | 58 6d 69 6e 3b 22 20 20 |"Xmin=";|Xmin;" |
|00002220| 20 22 3b 22 58 6d 61 78 | 3d 22 3b 58 6d 61 78 0d | ";"Xmax|=";Xmax.|
|00002230| 11 9e 2f ec 20 38 30 2c | 39 31 35 20 3a 20 f1 22 |../. 80,|915 : ."|
|00002240| 59 6d 69 6e 3d 22 3b 59 | 6d 69 6e 3b 22 20 20 20 |Ymin=";Y|min;" |
|00002250| 22 3b 22 59 6d 61 78 3d | 22 3b 59 6d 61 78 0d 11 |";"Ymax=|";Ymax..|
|00002260| a8 0f ef 20 34 20 3a 20 | e6 20 34 2c 33 0d 11 b2 |... 4 : |. 4,3...|
|00002270| 09 65 72 72 3d 30 0d 11 | bc 05 20 0d 11 c6 05 e1 |.err=0..|.. .....|
|00002280| 0d 11 d0 05 20 0d 11 da | 04 0d 11 e4 2f f4 20 2a |.... ...|..../. *|
|00002290| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000022a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000022b0| 2a 2a 2a 2a 2a 2a 2a 2a | 0d 11 ee 04 0d 11 f8 30 |********|.......0|
|000022c0| f4 20 55 73 65 72 20 73 | 75 70 70 6c 69 65 64 20 |. User s|upplied |
|000022d0| 72 6f 75 74 69 6e 65 20 | 74 6f 20 65 76 61 6c 75 |routine |to evalu|
|000022e0| 61 74 65 20 66 75 6e 63 | 74 69 6f 6e 0d 12 02 2b |ate func|tion...+|
|000022f0| f4 20 49 6e 70 75 74 20 | 20 70 61 72 61 6d 65 74 |. Input | paramet|
|00002300| 65 72 73 20 20 5a 31 3d | 52 65 7b 5a 7d 20 20 5a |ers Z1=|Re{Z} Z|
|00002310| 32 3d 49 6d 7b 5a 7d 0d | 12 0c 31 f4 20 4f 75 74 |2=Im{Z}.|..1. Out|
|00002320| 70 75 74 20 70 61 72 61 | 6d 65 74 65 72 73 20 20 |put para|meters |
|00002330| 55 31 3d 52 65 7b 46 75 | 6e 63 7d 20 20 55 32 3d |U1=Re{Fu|nc} U2=|
|00002340| 49 6d 7b 46 75 6e 63 7d | 0d 12 16 3c f4 20 53 65 |Im{Func}|...<. Se|
|00002350| 74 20 65 72 72 6f 72 20 | 70 61 72 61 6d 65 74 65 |t error |paramete|
|00002360| 72 20 65 72 72 20 74 6f | 20 30 20 66 6f 72 20 73 |r err to| 0 for s|
|00002370| 75 63 63 65 73 73 66 75 | 6c 20 65 76 61 6c 75 61 |uccessfu|l evalua|
|00002380| 74 69 6f 6e 0d 12 20 26 | f4 20 53 65 74 20 65 72 |tion.. &|. Set er|
|00002390| 72 3d 31 20 66 6f 72 20 | 69 6c 6c 65 67 61 6c 20 |r=1 for |illegal |
|000023a0| 61 72 67 75 6d 65 6e 74 | 20 5a 0d 12 2a 25 f4 20 |argument| Z..*%. |
|000023b0| 53 65 74 20 65 72 72 3d | 32 20 66 6f 72 20 66 75 |Set err=|2 for fu|
|000023c0| 6e 63 74 69 6f 6e 20 6f | 76 65 72 66 6c 6f 77 0d |nction o|verflow.|
|000023d0| 12 34 04 0d 12 3e 0b dd | 20 f2 66 75 6e 63 0d 12 |.4...>..| .func..|
|000023e0| 48 05 20 0d 12 52 09 65 | 72 72 3d 30 0d 12 5c 05 |H. ..R.e|rr=0..\.|
|000023f0| 20 0d 12 66 1c f4 20 20 | 3c 20 62 6f 64 79 20 6f | ..f.. |< body o|
|00002400| 66 20 70 72 6f 63 65 64 | 75 72 65 20 3e 0d 12 70 |f proced|ure >..p|
|00002410| 05 20 0d 12 7a 05 e1 0d | 12 84 05 20 0d 12 8e 32 |. ..z...|... ...2|
|00002420| f4 20 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |. ******|********|
|00002430| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002440| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 0d 12 |********|******..|
|00002450| 98 05 20 0d 12 a2 22 f4 | 20 45 78 61 6d 70 6c 65 |.. ...".| Example|
|00002460| 20 64 61 74 61 20 46 28 | 5a 29 20 3d 20 5a 5e 33 | data F(|Z) = Z^3|
|00002470| 30 20 2d 20 31 0d 12 ac | 04 0d 12 b6 37 f4 20 46 |0 - 1...|....7. F|
|00002480| 69 72 73 74 20 70 61 72 | 61 6d 65 74 65 72 20 6d |irst par|ameter m|
|00002490| 61 79 20 62 65 20 50 4f | 4c 59 20 28 66 6f 72 20 |ay be PO|LY (for |
|000024a0| 72 65 61 6c 20 70 6f 6c | 79 6e 6f 6d 69 61 6c 29 |real pol|ynomial)|
|000024b0| 0d 12 c0 3e f4 20 20 20 | 20 20 20 20 20 20 20 20 |...>. | |
|000024c0| 20 20 20 20 20 20 20 20 | 20 20 6f 72 20 46 55 4e | | or FUN|
|000024d0| 43 20 28 66 6f 72 20 75 | 73 65 72 20 73 75 70 70 |C (for u|ser supp|
|000024e0| 6c 69 65 64 20 66 75 6e | 63 74 69 6f 6e 29 0d 12 |lied fun|ction)..|
|000024f0| ca 04 0d 12 d4 38 f4 20 | 53 65 63 6f 6e 64 20 70 |.....8. |Second p|
|00002500| 61 72 61 6d 65 74 65 72 | 20 69 73 20 4e 72 20 74 |arameter| is Nr t|
|00002510| 68 65 20 6e 75 6d 62 65 | 72 20 6f 66 20 7a 65 72 |he numbe|r of zer|
|00002520| 6f 65 73 20 73 6f 75 67 | 68 74 0d 12 de 43 f4 20 |oes soug|ht...C. |
|00002530| 28 46 6f 72 20 61 20 70 | 6f 6c 79 6e 6f 6d 69 61 |(For a p|olynomia|
|00002540| 6c 20 74 68 69 73 20 69 | 73 20 6f 6e 65 20 6c 65 |l this i|s one le|
|00002550| 73 73 20 74 68 61 6e 20 | 74 68 65 20 6e 75 6d 62 |ss than |the numb|
|00002560| 65 72 20 6f 66 20 63 6f | 65 66 66 73 29 0d 12 e8 |er of co|effs)...|
|00002570| 04 0d 12 f2 40 f4 20 54 | 68 69 72 64 20 70 61 72 |....@. T|hird par|
|00002580| 61 6d 65 74 65 72 20 6d | 61 79 20 62 65 20 52 45 |ameter m|ay be RE|
|00002590| 41 4c 20 20 20 20 28 69 | 66 20 6f 6e 6c 79 20 72 |AL (i|f only r|
|000025a0| 65 61 6c 20 72 6f 6f 74 | 73 20 73 6f 75 67 68 74 |eal root|s sought|
|000025b0| 29 0d 12 fc 35 f4 20 20 | 20 20 20 20 20 20 20 20 |)...5. | |
|000025c0| 20 20 20 20 20 20 20 20 | 20 20 20 6f 72 20 43 4f | | or CO|
|000025d0| 4d 50 4c 45 58 20 28 74 | 68 65 20 75 73 75 61 6c |MPLEX (t|he usual|
|000025e0| 20 63 61 73 65 29 0d 13 | 06 04 0d 13 10 38 f4 20 | case)..|.....8. |
|000025f0| 54 68 65 20 6c 61 73 74 | 20 66 6f 75 72 20 70 61 |The last| four pa|
|00002600| 72 61 6d 65 74 65 72 73 | 20 58 6d 69 6e 2c 20 58 |rameters| Xmin, X|
|00002610| 6d 61 78 2c 20 59 6d 69 | 6e 20 61 6e 64 20 59 6d |max, Ymi|n and Ym|
|00002620| 61 78 0d 13 1a 37 f4 20 | 64 65 66 69 6e 65 20 74 |ax...7. |define t|
|00002630| 68 65 20 67 72 61 70 68 | 69 63 73 20 61 72 65 61 |he graph|ics area|
|00002640| 20 64 69 73 70 6c 61 79 | 65 64 20 62 79 20 74 68 | display|ed by th|
|00002650| 65 20 70 72 6f 67 72 61 | 6d 0d 13 24 04 0d 13 2e |e progra|m..$....|
|00002660| 3a f4 20 49 66 20 66 69 | 72 73 74 20 70 61 72 61 |:. If fi|rst para|
|00002670| 6d 65 74 65 72 20 69 73 | 20 50 4f 4c 59 2c 20 74 |meter is| POLY, t|
|00002680| 68 65 6e 20 74 68 65 20 | 6e 65 78 74 20 64 61 74 |hen the |next dat|
|00002690| 61 20 6c 69 6e 65 73 0d | 13 38 3a f4 20 63 6f 6e |a lines.|.8:. con|
|000026a0| 74 61 69 6e 20 74 68 65 | 20 4e 72 2b 31 20 63 6f |tain the| Nr+1 co|
|000026b0| 65 66 66 69 63 69 65 6e | 74 73 20 43 28 4e 72 29 |efficien|ts C(Nr)|
|000026c0| 2c 20 43 28 4e 72 2d 31 | 29 2e 2e 2e 2e 43 28 30 |, C(Nr-1|)....C(0|
|000026d0| 29 0d 13 42 04 0d 13 4c | 15 dc 20 50 4f 4c 59 2c |)..B...L|.. POLY,|
|000026e0| 33 30 2c 43 4f 4d 50 4c | 45 58 0d 13 56 17 dc 20 |30,COMPL|EX..V.. |
|000026f0| 2d 31 2e 32 2c 31 2e 32 | 2c 2d 31 2e 32 2c 31 2e |-1.2,1.2|,-1.2,1.|
|00002700| 32 0d 13 60 07 dc 20 31 | 0d 13 6a 19 dc 20 30 2c |2..`.. 1|..j.. 0,|
|00002710| 30 2c 30 2c 30 2c 30 2c | 30 2c 30 2c 30 2c 30 2c |0,0,0,0,|0,0,0,0,|
|00002720| 30 0d 13 74 19 dc 20 30 | 2c 30 2c 30 2c 30 2c 30 |0..t.. 0|,0,0,0,0|
|00002730| 2c 30 2c 30 2c 30 2c 30 | 2c 30 0d 13 7e 1a dc 20 |,0,0,0,0|,0..~.. |
|00002740| 30 2c 30 2c 30 2c 30 2c | 30 2c 30 2c 30 2c 30 2c |0,0,0,0,|0,0,0,0,|
|00002750| 30 2c 2d 31 0d ff | |0,-1.. | |
+--------+-------------------------+-------------------------+--------+--------+
